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Abstract 

We examine the bound state (s) associated with a single cubic nonlinear im- 
purity, in a one-dimensional tight-binding lattice, where hopping to first- 
and-second nearest neighbors is allowed. The model is solved in closed form 
via the use of the appropriate lattice Green function and a phase diagram is 
obtained showing the number of bound states as a function of nonlinearity 
strength and the ratio of second to first nearest-neighbor hopping parame- 
ters. Surprisingly, a finite amount of hopping to second nearest neighbors 
helps the formation of a bound state at smaller (even vanishingly small) non- 
linearity values. As a consequence, the selftrapping transition can also be 
tuned to occur at relatively small nonlinearity strength, by this increase in 
the lattice dispersion. 
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The effects of impurities on the transport properties of materials contin- 
ues to be an interesting subject. When the concentration of impurities inside 
a material is finite, we speak of a disordered system. In one dimensional 
systems, this disorder gives rise to the well-known phenomenon of Anderson 
localization, where all states acquire a finite localization length. This pre- 
cludes any amount of transport in the system. Economou and coworkers [l]] 
have argued that there exist a correspondence between the property that, in 
dimensions smaller or equal to two, all states of a disordered system are local- 
ized, and the property that an isolated impurity can always trap a particle. 
More recently, the problem of nonlinear impurities have received considerable 
attention. In a condensed matter context, they appear in strongly coupled 
electron-vibration systems, when the vibrational degrees of freedom have 
the ability of adapting quickly to the presence of the electron, giving rise to 
polaronic behavior 0. Nonlinear impurities appear also in other fields, such 
as nonlinear optics. For instance, an array of linear waveguides containing 
a single or several nonlinear, Kerr-like guides. The transversal dynamics de- 
scribing the energy exchange among waveguides is formally identical to the 
dynamics of an excitation propagating in a linear tight-binding lattice in the 
presence of one or more nonlinear impurities p3]. 

In systems where an electron (or excitation) is propagating while strongly 
interacting with vibrational degrees of freedom, an "effective" nonlinear evo- 
lution equation for the electron can be obtained, under the assumption that 
the vibrations adapt instantly to the presence of the electron^]. This equa- 
tion, known as the Discrete Nonlinear Schrodinger (DNLS) equation, has the 
form: 



where C n is the probability amplitude of finding the electron on lattice site 
n at time t, V is the nearest-neighbor hopping parameter and Xn is the 
nonlinearity parameter at site n proportional to the square of the electron- 
vibration coupling. In the conventional DNLS equation, the sum in (|]) is 
restricted to nearest-neighbors (n.n.). 

For the case of a single nonlinear impurity in a onedimensional lattice 
(Xn — ^n,o X)i we have obtained in a previous work|4| its bound state an- 
alytically, via lattice Green functions, and have shown that a bound state 





(1) 
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is possible provided |x/2V| > 1. This result deviated markedly from the 
well-known linear impurity case where, a bound state exists for any impurity 
strength. The extension to an impurity of arbitrary nonlinear exponent a 
(i.e., |C n | Q instead of |C n | 2 in ([!])), revealed^], ^| that for a < 2 there is 
always a bound state for any finite xjV . At a = 2 (the standard DNLS 
case) there is one bound state for \x/2V\ > 1, while for a > 2 there is a 
critical curve in x~ oc space, below which there is no bound state, while above 
it, there are two bound states. On the critical curve, there is a single bound 
state. 

When this nonlinear impurity is embedded in a square lattice^, the x~ 
a bound state phase diagram shows a single curve separating two regimes. 
Below the curve, there are no bound states; on the curve there is a single 
bound state, while above the curve there are two bound states. One of 
these become more localized upon increment of the nonlinearity parameter 
while the other becomes more delocalized. Bound states for single nonlinear 
impurities have also been computed for other systems, including a Cayley 
tree||, a triangular lattice ]10| and a cubic lattice || |10 |. 



In all the studies above, only dispersion to first nearest-neighbors has 
been considered. For systems with long-range dispersion, a continuum- like 
approximation that employs a nonlocal nonlinear Schrodinger equation has 
been proposed [TT|. For the case of a discrete system with a hopping param- 



eter of the form V nm = V/\n — m\ s , another study|12fl followed a variational 
approach, based on a plausible ansatz for the localized state. Among other 
things, they found that there is a critical s cr such that all dispersive in- 
teractions decreasing faster than r~ scr lead to similar qualitative behavior 
as the DNLS with only nearest neighbor transfer. In this work we solve 
in closed form via Green functions, the case of a single DNLS impurity in a 
tight-binding lattice, including hopping to first and second nearest -neighbors. 
Since the phenomenon of selftrapping is the result of the quantum struggle 
between the tendency to spread (dispersion) and the tendency to localize 
(nonlinearity), one might surmise that any increase in dispersion will have 
the simple effect of increasing the nonlinearity needed to selftrap. However 
as we will see, this is not necessarily the case and a small increment in disper- 
sion can actually favor the formation of a bound state at smaller nonlinearity 
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strength. 



1 Onedimensional lattice with dispersion 

Let us consider the problem of determining the existence of bound states and 
dynamic selftrapping characteristics for an electron (or an excitation) moving 
on a onedimensional dispersive lattice with hopping up to second nearest- 
neighbors, which contains a single DNLS impurity at the origin n — 0. The 
DNLS equation (|l|) reduces to 

i (^j = ^i(C n+1 + C n _ x ) + V 2 (C n+2 + C„_ 2 ) - x|C | 2 C 8 n>0 . (2) 
For stationary states, one puts C n (t) = exp(-iEt) n , obtaining 

E <p n = V 1 (4> n+ i + 0„_i) + V 2 (<p n+2 + (pn-2) - X \4>o\ 2 4>o <W (3) 
The Hamiltonian that gives rise to (0) is 

H = H + H U (4) 

where 

H = V 1 XX|n)(n+l| + |n + l)(n| ) + V 2 J2(\n)(n + 2| + \n + 2){n\ ) (5) 

n n 

and 

^i = -xl0o| 2 |O)(O|. (6) 

The {\n)} represent Wannier electronic states and ^(X^) is the nearest (next- 
to-nearest) neighbor transfer matrix element. In the absence of impurity, the 
energy band is given by 

E(k) = 2V 1 cos(k) + 2V 2 cos(2A;). (7) 

A simple analysis shows that, for positive V 2 /Vi, the upper and lower band 
edges obey 

E max = 2{Vi + V 2 ) (8) 

f _ 2 (^ - V 2 ) V./V, < 1/4 

a™ - j -(y?/W 2 ) - 2V 2 V 2 /V! > 1/4 [J) 
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As a result, while the upper edge always increase linearly with V2, the 
lower edge first decreases (in magnitude), reaching a minimum value of 
2(1 - (l/v^))^ « 1.29Vi at V 2 = (l/y/S)Vi. Afterwards, the lower edge 
increases in magnitude with V 2 . At large V%, this increase will be almost 
linear. These features will of importance in the next section, where we de- 
termine the position of the impurity bound state (s). 

The formalism of lattice Green functions is particularly suitable for the 
problem of determining the bound state (s). First, we normalize all energies 
by V x and define: z = E/V x , H = H/V u 7 = x/Vi and 5 = V 2 /V 1 . The di- 
mensionless lattice Green function SG = 1/ (z—H) can be formally expanded 
as[]l4 

G = G<® + G®H X G® + G^HxG^H^ + ... (10) 

where G^ is the unperturbed (7 = 0) Green function and Hi = — 7|0o| 2 |0 >< 0| 
The sum in Eq.(^) can be carried out exactly to yield 

r — r>(°) m0 0n (-\-\\ 



l + l\M 2 G 



oo 



where G mn = (m\G\n) and 

G (0) (z) - J_ f U exp^(m - n) ] 

whose evaluation in closed form is done in the Appendix. 
1.1 Bound States 



The energy of the bound state(s), Zb is obtained from the poles of G mn , i.e., 
by solving 1 = — 7)00^ 1 2 Coo • The bound state amplitudes <fi$ are obtained 
from the residues of G mn (z) at z — z b . In particular, 



#| 2 = Res{GooW}_ b = -^yf^ (13) 

^00 \ z b) 



Inserting this in the bound state energy equation leads to 

}_ _ ^00 ( z b) 

7 G q (z b ) 



(14) 
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Now, as the reader can easily verify using (|HD and (|T2|), changing both, the 
sign of xjV\ and the sign of V2/V1 has the effect of changing the sign of the 
corresponding eigenvalue. Thus, for a complete parameter space examination 
of all the possible eigenvalues, we only need to consider a fixed sign for V 2 /Vi 
(say, positive) and the two possible signs of xjV\. 

Figure 1 shows the left- and right-hand side of (14) for several (positive) 
values of 5. When 7 > (Fig. la) and as S increases from zero, the RHS 
of (0) moves towards the origin and increases its height, until S reaches 
1/4, where the height diverges. Further increase in S decreases the height 
of the curves, but they continue to approach the origin until 5 = 1/ y/8. 
Afterwards, the curves move away from the origin while their heights continue 
to decrease. For the case of 7 < (Fig. lb), the situation is quite different: 
For a given 5 value, there is a minimum I7J at which there is a bound state. 
Further increase in I7I creates two bound states, one of which will ultimately 
disappear at a further finite |7&| value, leaving only a single bound state. 

Figure 2 displays a phase diagram in nonlinearity-dispersion space show- 
ing the number of bound state(s). For positive nonlinearity, the critical curve 
separating the region with no bound states from the region with one bound 
state, decreases with increasing 5 and reaches zero at 5 = 1/4. Afterwards, 
it increases monotonically with further dispersion. Thus, there is a finite 
dispersion interval, < 5 < 1/4 where, contrary to what might be expected, 
an increase in dispersion actually reduces the nonlinearity needed to create a 
bound state. This can also be seen in Fig. 3a, which shows the bound state 
energy as a function of (positive) nonlinearity, for several values of dispersion 
5. This reduction in nonlinearity needed to sustain a bound state is, of course 
due to the reduction in the width of the negative portion of the band with 
a small positive dispersion, and thus, it has a linear origin. In the negative 
nonlinearity sector, we have in Fig.2 two critical curves separating regions 
with no bound states, two bound states and one bound state. Here, an 
increase in dispersion, causes a corresponding increase in the minimum non- 
linearity needed for the creation of a bound state(s). This is correlated with 
the fact that the width of the positive portion of the band always increases 
linearly with 5. Figure 3b, shows the bound state energy as a function of 
nonlinearity, for negative 7. Here, for a given value of S, there exists a critical 
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nonlinearity value 7 a for which there is a bound state, with energy outside 
the band. Further increase in nonlinearity creates two bound states, one of 
which increases its energy monotonically with nonlinearity while the other 
state decreases its energy towards the band, reaching it at a finite nonlinear- 
ity value 7b. As Fig. 4 shows, in the regime with two bound states, as the 
magnitude of the nonlinearity is increased, one of the states becomes more 
localized on the impurity site, while the other becomes more delocalized, 
ultimately disappearing into the continuum at a finite nonlinearity value. 

1.2 Transmission across the impurity 

Inside the band, z m i n (8) < z < z max {5) (see (||) and (0)), all states are 
extended and given by| 



m = \k) + G^ + T + (z) \k), (15) 
where |A;) is a plane wave and the scattering-matrix T is 

T = H 1 + H 1 G {0) H x + #i G (0) E x G (0) H x + ... (16) 
This can be summed exactly, to yield 

T= ^mm. (17) 

i + tICoI 2 g£> 

With this, we can compute the scattering amplitude at any site n. In par- 
ticular, the scattering amplitude at the impurity site Co = is 

c ,. M -«™«. (18) 

-L + Tl^ol ^oo 

The transmission coefficient is then given by the probability at the impurity 
site, t = \C \ 2 . From (Tl8|), we have 



|l+ 7 f Gffil' l + 7 2 * 2 Im[Gto] 



(19) 



which leads to the cubic equation: 7 2 t 3 Im[Goo(<2)] 2 + t — 1 = 0. This i 



is 



invariant under 7 — > —7, implying that the transmission does not depend on 
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the sign of the nonlinearity parameter. The physical solution for t is 

2 

-2 6^ + (l8 7 Im[Gto] + 2 V 4 + 27 7 2 Im[G J (z)] 2 ) * 
i = v / _^ ^20) 

2 3i 7 MGS(*)] (97MGSW] + v / 3 V /4 + 277 2 Im[GS(^)] 2 )" 

The specific form of ImfGgo (2)] in our case, is given in closed form in the 
Appendix. Figure 5 shows several transmission curves as functions of the 
plane waves dimensionless energy z, for several different ratios 5. The most 
remarkable new feature is the appearance of an abrupt "dip" on the trans- 
mission near the lower edge of the band at 5 ~ 0.4. As 5 increases further, 
the "dip" moves to the right and eventually (not shown) approaches the up- 
per band edge and merges with it. This "dip" is related to the creation of a 
secondary "branch" in Im[Goo (z)] ; as shown in Fig. B of the Appendix. 

1.3 Dynamic Selftrapping 

We place the electron at the impurity site at t = and observe its time 
evolution, according to Eq.fl2|). The observable of interest here is the long- 
time average probability of finding the electron on the initial site after a 
relatively long time T: 

P = lim(l/T) f T \C (t)\ 2 dt, |C o (0)| = l. (21) 

By using the transformation C n — > (— l) n C n , that leaves \C n \ 2 invariant, 
the reader can verify that changing V2 — * —V2 and x ~^ ~Xi transforms 
(0) into its complex conjugate. Thus, for our completely localized real initial 
condition Co(0) = 1, the observable Pq is invariant under the transformation. 
Therefore, for a complete parameter study, it is enough to consider a positive 
V2 (or S) and the two possible signs for the anharmonicity \ ( or 7)- 

We use a fourth-order Runge-Kutta numerical scheme, whose accuracy 
is monitored through total probability conservation: 1 = J2n\C n (t)\ 2 - To 
avoid undesired boundary effects, a self-expanding lattice is usedfl3|]. Figure 
6 shows Po as a function of (positive) nonlinearity parameter 7, for several 
different dispersion 5 values. As anticipated from the bound state results, 
an increase of dispersion reduces the critical nonlinearity for the onset of 
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selftrapping. The minimum threshold occurs around 5 ~ 0.3. In the im- 
mediate vicinity of this value, the selftrapping transition also seems to lose 
some sharpness. Subsequent increase in 5 increases the critical nonlinearity 
7 C again and restores sharpness to the Pq curve. At 5 ~ 1, the P curve 
almost coincides with the 5 = case. Thereafter, 7 C increases in an almost 
linear fashion with 5. For the case of a negative nonlinearity parameter (not 
shown), the critical nonlinearity always increases monotonically with disper- 
sion. As explained above, this case corresponds also to 7 > 0, 5 < 0. 



2 Discussion 



We have examined the conditions for the formation of a bound state at a 
nonlinear impurity site, in a onedimensional linear lattice with hopping to 
first and second nearest-neighbors. The formalism employed lattice Green 
functions, which have been evaluated in closed form for our system. We 
found that this increment in dispersion can actually favor the formation of 
a bound state at smaller nonlinearity strength. As a consequence, the onset 
of dynamical selftrapping at the impurity site can also be shifted to lower 
nonlinearity thresholds. This tuning effect could have some impact on the 
design of completely nonlinear optical or solid nanostructures, where one is 
interested in static ( "selftrapped" ) or mobile discrete nonlinear excitations 
("discrete solitons"). In the former, one is interested in achieving abrupt 
selftrapping at relatively low nonlinearity (or input power) parameter values; 
in the latter, the probability profile corresponding to an impurity bound 
state, could be useful as an initial condition for the launch of a discrete 
soliton, which can carry energy across a nonlinear fiber array from a given 



guide to another, distant one [15 



3 Acknowledgments 

This work was supported in part by FONDECYT grant 1020139. 



9 



4 APPENDIX: THE LATTICE GREEN FUNC- 
TION 



The dimensionless lattice Green function for our tight-binding chain including 
dispersion up to second nearest-neighbors is 

r(°) (?) - 1 f M exp[z0(m - n) ] 

Umn[Z) ~ 2vr L, d<P [z- cos(0) - 5 cos(20) ] • [ZZ) 



We follow the usual method [|T| of transforming fl2"2] ) to an integral over a 
complex variable w along the unit circle. Thus, we have 

G (o) (z) = _J_ / *» W "' m ' +1 (23) 

mnl J 2ttz/ 5w 4 + - w 2 2 + u> + 5 1 J 

The roots of 5 u> 4 + w 3 — u; 2 ,2 + u> + 5 = are given by 



Z-2 



2.3 



2.4 



-1 + 5 71 ( 


» 


-A 


'7a( 


» 


-5 7iW) 






45 








-l-5 7 i( 


» 


+ A 


^72 ( 


» 


+ 5 7 i(z)) 






45 








-(1 + ^71 




) + A 


/72 




+ *7i(*)) 






45 








-1+5 7l ( 


» 


+ A 


/7a ( 


» 


-5 7^)) 


45 



where 71(2) = ^8 + (1/5 2 ) + (4z/5) and 72(2) = 1 - 45 2 + 25;z. Only the 
poles contained inside the unit circle contribute to the integral. Thus, 



\zi\ < 1 =>• 2(1 + 5) <z< 00 

[ -00 < z < 00, 5 < 1/4 

\ -00 < z < -1/(45) - 25 and 2(5 - 1) < z < 00, 5 > 1/4 

I #3 1 > 1 =>- no contribution 

J -00 < 2 < -2(1-5), 5 < 1/4 
) -00 < z < -1/(45) - 25, 5 > 1/4 



N < 1 



IzA < 1 



Using the above, we can express G^(z) as 



27r5 J 5(u> — Zi)(u> — z 2 )(w — z 3 )(w — 24) 
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e[*-(2(i + *)] 4'" m|+1 



<5(^i - 22X21 - ^)(^i - ^4) 

[6((l/4) - 5) + 0(5 - (l/4))[0(-l/(45 - 25 - 2z)) + 9(z - 2(5 - 1))]] zjf""'* 1 



5(z 2 - z 1 )(z 2 - z 3 )(z 2 - z 4 ) 
[0((l/4) - 5)0(-2(l -S)-z) + 0(5 - (l/4))0(-l/(45) - 25 - 2z)\ 4" m|+1 



+ 



5(^4 - ^l)(^4 - ^2) (^4 - Z 3 ) 

After some tedious algebra, we find for the real and imaginary parts of the 
diagonal Green function: 



(25) 



Re[G ( S(z 



where, 



and 



' 2V25 [(7i5 + 12) Jl - 5 (45 - 2z + 71)+ 



(7i* " 72) Jl + 5(-45 + 2z + 7l ) J- 1 z > ^(5) 



- 2V25 [( 7l 5 + 72)^1-5 (45-2z + 7l )- 

(7i* - 72) + 6 (-45 + 2z + 7l ) ]- x ^<z min 



otherwise 



(26) 



-2(1-5) 5<l/4 
-1/(45) -25 5 >l/4 

na:E = 2(1 + 5) 



Jro[Gffi(*)] = 



0((l/4) - 5)[G(z - 2(5 - 1)) - Qjz - 2(1 + 5))] | 

2^1 -di(z) 2 |l + 45 
0(5 - (l/4))[6(z + (1/(45)) + 25) - G(z - 2(1 + 5))] 



+ 



2yjl -d^z) 2 ] 1 + 45 di^)! 
0(5 - (l/4))[8(z + (1/(45)) + 25)) - 8(* - 2(5 - 1))] 
2 v /l-4W 2 |l + 45 d 2 (*)| 



(27) 



where, ^(2) = (1/4 5)(-l + ^/l + 8 5 (5 + ^/2))) and d 2 (z) = (1/4 5)(-l- 
^1 + 8 5(5 + ^/2))). 

Figures A and B, show Re[G^(z)\ and lm[G^Q (z)], respectively for different 
values of dispersion 5. 
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Captions List 



Fig.l : Left-and-right hand side of the eigenenergy equation ([14]) , for several 
dispersion 5 values, for positive (a) and negative (b) nonlinearity. 

Fig. 2 : Bound state phase diagram in nonlinearity-dispersion space (7 = 
X /V 1 ,5 = V 2 /V 1 ) 

Fig. 3 : Bound state energy as a function of positive (a) and negative(b) 
nonlinearity, for several dispersion values. 

Fig. 4 : Probability at the impurity site as a function of nonlinearity, in the 
negative nonlinearity sector, for several dispersion values. 

Fig. 5: Transmission coefficient of plane waves across the nonlinear impurity 
as a function of the plane wave dimensionless energy, for several dispersion 
values. 

Fig. 6: Long-time average probability of finding the electron on the impurity 
site, as a function of nonlinearity, for several different dispersion values. (T = 
203 Vi) 

Fig. A: Real part of the diagonal unperturbed Green function Gqq(z), for 
several dispersion values: (a), 0.2(b), 0.4(c), 0.6 (d), 0.8 (e) and 1.0(e). 

Fig. B: Same as in Fig. A, but for the imaginary part of Gqq(z). 
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